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A momentum balance equation is developed to investigate the magnetotransport properties in 
monolayer molybdenum disulphide when a strong perpendicular magnetic field and a weak in¬ 
plane electric field are applied simultaneously. At low temperature, in the presence of intravalley 
impurity scattering Shubnikov de Haas oscillation shows up accompanying by a beating pattern 
arising from large spin splitting and its period may halve due to high-order oscillating term at large 
magnetic field for samples with ultrahigh mobility. In the case of intervalley disorders, there exists a 
magnetic-field range where the magnetoresistivity almost vanishes. For low-mobility layer, a phase- 
inversion of oscillating peaks is acquired in accordance with recent experiment. At high temperature 
when Shubnikov de Haas oscillation is suppressed, the magnetophonon resonances induced by both 
optical phonons (mainly due to homopolar and Frohlich modes) and acoustic phonons (mainly due 
to intravalley transverse and longitudinal acoustic modes) emerge for suspended system with high 
mobility. For the single layer on a substrate, another resonance due to surface optical phonons 
may occur, resulting in a complex behavior of the total magnetoresistance. The beating pattern of 
magnetophonon resonance due to optical phonons can also be observed. However, for nonsuspended 
layer with low mobility, the magnetoresistance oscillation almost disappears and the resistivity 
increases with field monotonously. 

PACS numbers: 75.47.-m, 72.20.-i, 81.05.Hd 


I. INTRODUCTION 

The rise of graphene^!*! showing outstanding mechan¬ 
ical and electronic properties, launched the era of mono- 
layer material. However, pristine graphene does not have 
a band gap, a property essential for electronic applica¬ 
tions. Although, it is possible to open small band gap 
in graphene by some method)^ it will inevitably lead 
to increased fabrication complexity and reduced perfor¬ 
mance of devices££ This produces great limitations on 
its becoming a perfect candidate for the next genera¬ 
tion nanoelectronic material. In contrast to graphene, 
the transition metal dichalcogenides are semiconductors 
with naturally occurring band gap, which overcomes this 
problem directly. A prominent representative in this 
dichalcogenide family is molybdenum disulphide (M0S2). 
Bulk M0S2 has an indirect gap, while monolayer M0S2, 
which can be isolated by exfoliation techniques similar 
to graphene, is a direct-gap semiconductor with a gap 
of 1.9 eVi! Due to the large carrier mobility!, high cur¬ 
rent carrying capacity!, strong spin-orbit coupling, and 
coupling of spin and valley degrees of freedom, mono- 
layer M0S2 may become a replacement of graphene or 
even a candidate for the exploitation of novel valleytronic 
devices »!! 

On the aspect of transport investigation of monolayer 
M0S2, the linear mobility is close to 200cm 2 /Vs at low 
temperature where a high-re gate dielectric was used to 
suppress the charged-impurity scattering strongly: 11 - This 
value is still lower than the theoretical prediction, where 
the highest phonon-limited mobility in n-type mono¬ 


layer M0S2 is 410 cm 2 /Vs at room temperature*!! On 
the other hand, the single layer M0S2 device grown by 
chemical vapor deposition shows low temperature mo¬ 
bility up to 500cm 2 /Vs, where the leading scattering 
mechanism is believed to be the short-range scatter- 
ers at high carrier density*!! Hence, the main scatter¬ 
ings determining the linear mobility is still open ques¬ 
tion. Further, looking at all theoretical studies on electric 
transport , 12 i 14 i 15 the strong spin-orbit coupling in n-type 
monolayer M0S2, which can lead to interesting coupled- 
spin-valley physics ) 10 ! 16 ’ 17 is omitted completely and the 
energy band is chosen to be a simple parabolic one. 

Especially, in magnetotransport the spin-orbit cou¬ 
pling is important, which may result in the beating pat¬ 
tern of Shubnikov de Haas oscillation (SdHO)!! and in¬ 
duce direct magnetoresistance oscillation*!! Due to the 
spin-valley coupling, the magnetic control of the val¬ 
ley degree of freedom in monolayer M0S2 in the pres¬ 
ence of normal magnetic field has been achieved*!! The 
magneto-optical properties!! and magnetocapacitances!! 
have been analyzed in this system. However, even the ba¬ 
sic SdHO considering all kinds of scattering mechanisms 
in this single layer has not been seriously involved either 
in theoretical or in experimental works. Only recently, 
the SdHO was observed experimentally for the first time 
in monolayer and few-layer M0S2*!! In this paper, we ap¬ 
ply a momentum balance equation to investigate the lin¬ 
ear magnetotransport at both low and high temperatures 
including SdHO and magnetophonon resonance (MPR) 
effect induced by optical and acoustic phonons for both 
suspended and nonsuspended samples. 
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II. BASIC FORMULATION 

We consider a monolayer of transition metal dichalco- 
genide M 0 S 2 having large number N carriers in the x-y 
plane. These carriers, in addition to interacting with each 
other, are scattered by the random impurities and cou¬ 
pled with phonons in M 0 S 2 and substrate. There exists 
an external magnetic field B = (0, 0, B) applied along the 
z direction and a uniform electric field E in the mono- 
layer plane. The total Hamiltonian of the system is given 
by 


H e = H cm + Her, with 

Hem —"T7 ^ ' at(TjII x ajx T ByCjy') -f- NeE • R , 


=V n + NeE H, 


Her — yi 


. , , A 

%t[Tj 7 Tj x aj x + 7Tj y (Jj y J + jz 


XTj Sj z 


~ 1 


i<j 


(3) 

(4) 


H =H e +Hph+Hei+Hep. (1) 

Here the carrier part H e = J2j(^j+ er j' E )+Y^i<j Vc(r z - 
rj) with Tj = ( Xj,yj ) being the in-plane coordinate for 
the jth carrier and —e denoting its charge, and V c (ri — 
rj) standing for the Coulomb coupling potential between 
the ith and jth carriers, which, though depends solely 
on ri — rj for the thin-layer system, may vary with the 
spatial (^-direction) dielectric environment around the 
layer. The single-carrier low-energy Hamiltonian near K 
(-K) point in the Brillouin zone for the jth carrier in 
the presence of magnetic field is given b y 10 i 17 

hj = at{Tj'Kj x (jj x + njyCTjy ) + — <7j Z — XTjSj Z ® ^ . 

( 2 ) 

Here a is the lattice constant, t is the hopping integral, 
A is the energy gap, A is the spin-orbit coupling param¬ 
eter, Tj = ±1 is the valley index of the jth carrier re¬ 
ferring to ±K valley, 7 Tj = pj + eA(rj ) = (- nj x ,-Kj y ) 
is the canonical momentum with pj = ( Pj X ,Pj y ) be¬ 
ing the momentum of the jth carrier and the vec¬ 
tor potential in the Landau gauge A(rj) = {—Byj, 0), 
(Tj = (c Tj x ,aj y ,CTj Z ) is its pseudospin operator acting on 
the orbit {d z 2, (d x 2_ y 2 + iTjd xy )/V 2} and Sj Z is the 21- 
component of its real spin operator. There are Q val¬ 
leys locating near the halfway points along the F-K 
axes, which may introduce additional intervalley scatter¬ 
ing process, and thus influence the carrier transport^ 
However, the accurate value of the energy separation be¬ 
tween the K and Q points is still unsettle d 12 ’ 15 ! 24 with 
the estimate larger than 50meV. In the present calcu¬ 
lation, the Fermi energy of this discussed system is far 
below Q valleys. Hence, we can safely neglect the lim¬ 
ited effect of Q valleys and assume the system in the 
K- valley dominated carrier transport regime. H p h, H c i 
and H ep are phonon Hamiltonian, carrier-impurity and 
carrier-phonon interaction, whose forms can be found in 
the textbook^ and Rcfs.l26l and l27l. 

In terms of the center-of-mass (c.m.) momentum and 
coordinate defined as P = pj and R = N _1 Y^j r j 
for the whole system of N carriers and the relative-carrier 
momentum and coordinate p'j = pj — P/N and r' = 
Tj — R of the jth carrier t 27 i 28 the carrier Hamiltonian H e 
of this coupled many-body system can be written as the 
sum of a c.m. part H cm and a relative carrier part H er , 


Here IT = P+NeA(R) = (I7 Z , I7 y ) is the c.m. canonical 
momentum of the total system, n ' = p' + eA(r') = 
(7r' x , 7 Tj y ) is the canonical momentum for the jth relative 
carrier, and 

V = R = ~i[R, H] = al{T i a iB + Vjvh ( 5 ) 

j 

is the c.m. velocity operator of the carrier system. 

Note that, the commutation relation between the c.m. 
part Hem and the relative-carrier part H cr is of order 
of 1/N. Hence for a macroscopically large N system 
the c.m. motion and the relative motion of carriers are 
truly separated from each other. A spatially uniform 
electric field E shows up only in the c.m. part H cm , and 
Her is just the Hamiltonian of a monolayer M 0 S 2 sub¬ 
ject to a perpendicular magnetic field without the electric 
field. The coupling of two parts appears only through the 
carrier-impurity and carrier-phonon interactions. 

To proceed the calculation of transport properties in 
monolayer M 0 S 2 in the presence of a magnetic field, we 
can write down all the physical quantities in the Landau 
representation. The Landau levels of the single-particle 
Hamiltonian h is labeled by a band index a = ±1 for 
conduction and valence band, valley index r = ±1 for K 
and —K valley, and spin index s = ±1 for spin up and 
spin down in addition to the Landau index n with the 
form 


Earns = TsX + Ot\J (A - TsX)" + nU%, 
for n = 1, 2, 3 • • •, while for n = 0 

ErOs = -T (A - sX) + sX, 


( 6 ) 

(7) 


with A = A/2, A = A/2, and the cyclotron frequency 
uj c = \/2at/lB = \/2\e\Bat. One should take notice of 
the fact that the zero level ( n = 0) for K valley (r = 
+1) is in the valence band, while the zero level (n = 
0) for — K valley (t = —1) is in the conduction band. 
The corresponding eigenstates, including zero levels (n = 
0), are expressed as T arrls = \s ® k x ), with \s 

standing for the eigenstate of s z and 
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VnU 1 (r,k x ) = 


0 ik x x 


0r 


^rijfex (y) 

A n:rVn-i,fcx (y) 


(9) 


Here k x is the axcomponent of wave vector k, the coeffi¬ 
cient 


A“'I = 


\Jnu 


(A — rsA) — a.T\J (A — rsA) 2 + : 


( 10 ) 


and 0“’J = (A“’J) 2 + 1. Note that for K valley r = 
+1 {—K valley r = —1), only the valence band a = 
— 1 (conduction band a = +1) is allowed when n = 0. 
(f) n fc (y) is the harmonic oscillator eigenfunction giving 

by’ 


^n,fe x (y) 
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^2 n n\l B ^/n 


exp 


(y - yc ) 2 
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H n 
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(n) 

with H n (x) the Hermite polynomial, and y c = k x /(eB). 
In the Landau representation, the carrier-impurity and 
carrier-phonon Hamiltonians including both intravalley 
and intervalley interactions have the following forms: 


'Hoi =E E U TT ,(q)JZ' T lfs' {q y q . { R 


~r a ) 


Q’.T.n.S 


c arns C QL'r'n' s' 5 ( 1 ^) 

^ P =E E 


q,i/ a,r,n,s 
a'r'.n'.s' 


X C, 


arns^a'r'n's' • 


(13) 


Here U TT '(q) and M TT i(q, v) are the intravalley or inter¬ 
valley carrier-impurity scattering potential with r a being 
the impurity position and carrier-phonon coupling matrix 
of v branch, respectively; c aTns and c' aTns are the annihi¬ 
lation and creation operators of carrier; (j> qv = b qv + b_ qv 
is the phonon Held operator with b qu and b' qv being the 
annihilation and creation operators for a two-dimensional 
(2D) phonon of wave vector q in the branch v having fre¬ 
quency tt qv ; and the integral 


ja r n s 
J arns 


{q) = / dr' (<p%Ur',k x ) 


Vn'll' ( r \k x )). 


(14) 

The derivation of momentum balance equation starts 
from the rate of change of the c.m. canonical momentum 
n = -i[n,n\. To linear order in the carrier-impurity 
and carrier-phonon couplingsthe statistical average 
of this operator equation can be obtained by using the 
initial density matrix po = Z~ 1 e~^ n p h +' Her )/ T at temper¬ 
ature T in the case of weak in-plane electric field E. In 
the dc steady state, (II) = 0, the momentum balance 
equation for a system of unit area (N is thus understood 
as the carrier number density) reads 


0 = - Nev x B - NeE + / ei + f op , (15) 


with v = (V) being the averaged carrier drift velocity. 
The frictional forces experienced by the center of mass 
due to impurity and phonon scatterings, f OI and / ep , have 
the following form: 


/ei = n;^ \U TT '{q)\ 2 qn^ T ' (q,u 0 ), (16) 

q.r,r' 

/ ep = ^2 \M TTl (q,v)\ 2 qW 2 T \q, Q qv +w 0 ) 


^ H qv ^ ^ ^ f2 qu T u^o ^ 


(17) 


In the above expressions, rii is an effective impurity den¬ 
sity; n(x) = (e* — 1) _1 is the Bose distribution function; 
w 0 = q ■ v; HY (q,oj) is the imaginary part of the Fourier 
spectrum of the valley-dependent relative-carrier density 
correlation function, defined by 


H TT W - f) = -i0(t - f) ( [p T q T \t), pt T q (H)] } Q , (18) 

where p q T ( t ) = e* Wert Pq T e~ lHert with 

n TT ‘ = V T a ' T ' n ' s ' (a)Y r , , , , 

rq / j °arris v” ) L 'arns L '<x r'n's' 5 

a,n,s 
a ,n ,s 


and (• • • )o stands for the statistical averaging with re¬ 
spect to the initial density matrix p 0 i 27 i 28 

In most cases the electron density-correlation function 
in the presence of intercarrier coupling, IVY (q,ui), can 
be obtained in the random-phase approximation through 
the density-correlation function (q, u>) in the absence 
of intercarrier coupling, 


n Y'(q^) 


|e TT /(q,w)| 2 ’ 


(19) 


where e TT i(q,u) is the carrier-coupling related RPA 
screening function or carrier screening function, which 
may vary with the dielectric environment of two- 
dimensional (2D) monolayer. Therefore, in Eas. (11611 and 
(TlTl) n Y (<z,w) function can be replaced by njj (g,w) 
function, as long as the impurity and phonon scatter¬ 
ing potentials are considered screened by the intercarrier 
coupling: U TT '(q)/e TT >(q,u) and M TT >(q, v)/e TT >(q, w). 

The njj (q, w) function can be expressed as 


a 02 


(g,w) = 


27 
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with ImGo, rrls (e) standing for the imaginary part of re- while for intervalley case it has a more complex form 
tarded Green’s function G aTns (e) and the form factor for 
intravalley case is given by 
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with L™{z) being associated Laguerre polynomials, 0 = 
l^q 2 / 2, rii = min(n,n / ), 712 = max(n,n / ), mi = min(n — 
1 ,n'), m 2 = max(n — 1 ,n'), ki = min(n, n' — 1), fc 2 = 
max(n,n / — 1), 9 q is the polar angle of wave vector q , 
and 


f 1, n — 1 < n' 
y —1, n — 1 > n' 


( 1, n < n' — 1 
\ —1, n > n' — 1 


The momentum balance equation (1151) combining with 
equation of carrier density (l25l) describes the steady-state 
magnetotransport of monolayer M 0 S 2 , which can deter¬ 
mine either the drift velocity (charge current density) for 
given electric field or the electric field for given current. 
In the Hall configuration, e.g., with the charge current 
J (or drift velocity) in the x direction, J = (J, 0) = 
(— Nev,0), the momentum balance equation (1151) gives 
a transverse magnetoresistance R xy = —E y /(Nev) = 
—B/{Ne ) and a longitudinal magnetoresistance R xx = 

~(fei + fep)/{N 2 e 2 v). 


In the presence of carrier-impurity, carrier-phonon, 
and carrier-carrier scatterings, the Landau levels of 
monolayer M 0 S 2 are broadened. The imaginary part of 
the retarded Green’s function ImGo, Tns (e) or the density 
of state of the arnsfh Landau level is modeled using a 
Gaussian formal 


T ( \ 

Imb QTns (e) = -exp 

7 rvTr)..<? 


2(e Yarns') 


r 2 


(24) 


with r aTns denoting the half width. 

The chemical potential ey at temperature T is deter¬ 
mined by the carrier density (electron density N + or hole 
density 1V_) of the system by the following equation 


{ 


N + 



/(e)ImG +T „ s (e) 

[1 ^ /(e)] ImG_ Tns (e) 

(25) 


Here /(e) = {exp[(e — £f)/T] + l} -1 is the Fermi dis¬ 
tribution function. For electron conduction case, the 
summation index n in the above equation is taken over 
1, 2,3, • • • for K valley (r = +1), but 0,1, 2, • • • for —K 
valley (r = — 1). However, for hole conduction, it is taken 
over 0,1, 2, • • • for K valley, but 1, 2,3, • ■ • for —K valley. 


III. NUMERICAL RESULTS AND DISCUSSION 

For numerical calculation, we concentrate on the n- 
doped case, i.e., carrier is electron, N = N + and we only 
need to consider a = a 1 = +1. In the following, the index 
a or a' will be omitted. The half-width r rns should vary 
with the band indices generally. However, for simplicity, 
we neglect the effect of spin-orbit interaction, and take it 
with the form i 32 * 33 


r = 


eu} c oar 

7TTO* /r 


(26) 


Here /r is the zero-field mobility at temperature T, oo c 0 = 
eB/m* is the cyclotron frequency with effective mass 
m* = A/(a 2 i 2 ), and or is a phenomenological param¬ 
eter to relate the single-particle lifetime to the transport 
scattering time/ 32 ! 33 In the following numerical evalua¬ 
tion, we will set op = 3, except otherwise specified. 

The intravalley electron-impurity scattering potential 
is considered due to charged impurities distributed at a 
distance d from the layer 


U TT {q) 



2sonq 


(27) 
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with Zi standing for the effective impurity charge num¬ 
ber, k for the dielectric constant of MoS 2 . For suspended 
MoS 2 monolayer, d = 0. The carrier screening is taken 
into account with a static screening function of Thomas- 
Fermi for m 14 ! 31 ) 35 

eff 

e(q,0) = e(q) = l + ^. (28) 

For suspended M 0 S 2 , q x* equals the zero-temperature 
Thomas-Fermi wave vector qpp = m* e 2 / (tt£qk)', while 
for the layer on a substrate, the value of g^p, which could 
be a couple of times larger or smaller than qtf depending 
on the dielectric environment and carrier density, will be 
taken from Ref. [35]. Note that the main role of carrier 
screening is to enhance or decrease the mobility with or 
without magnetic field. In the present study, the effective 
impurity charge density riiZ 2 , which may be modified by 
the spatial dielectric environment of the system, is deter¬ 
mined by the zero-temperature carrier mobility qo in the 
absence of the magnetic field under the same screening 
condition, thus the major magnetic-field related behav¬ 
iors are not sensitive to the detailed form of the scattering 
potential or screening. 

In addition to above intravalley impurity scattering, we 
also include intervalley disorder scattering, which can be 
induced by lattice vacancy in a two-dimensional honey¬ 
comb lattice^ and by defects raised from ion irradiation 
as in graphene— The scattering potential is usually mod¬ 
eled by a d-function form, i.e., a constant U T f(q) = uq. 

For intrinsic electron-phonon couplings in suspended 
layer, we consider both intravalley and inter valley acous¬ 
tic deformation potential interactions. In the case of 
optical deformation potential, both zero-order and first- 
order couplings are taken into account and the homopo- 
lar mode is also included. The relevant formulas can be 
found in Ref. The polar longitudinal optical phonons 
are also important and their coupling matrix element 
with 2D carriers can be written a o 12 ’ 38 

M TT (q , Fr) = g F rerfc(gcr/2), (29) 



FIG. 1: (Color online) Longitudinal magnetoresistance as a 
function of average filling factor uo at different lattice tem¬ 
peratures T = 0.5, 1.0, 1.5, 2.5, 4.0 K. The inset shows the en¬ 
larged magnetoresistance at T = 0.5 K for small filling. The 
linear mobility at zero temperature^ q 0 = 4000 cm 2 /Vs and 
the electron density N = 7 x 10 12 cm~ 2 . 

f \ qv rather than the static function e(q, 0 ), should be 
used in the equation with phonon scattering. It has been 
show n 27 ! 42 that optic (as well as acoustic) phonon in¬ 
duced 2D resistivity with dynamic screening are essen¬ 
tially equivalent to those without screening at tempera¬ 
ture T > 100 K, when phonon scatterings play important 
roles. Therefore, in the numerical calculation of phonon- 
related magnetotransport at higher temperatures we will 
not include screening in the electron-phonon matrix ele¬ 
ment. 

The relevant parameters used in the numerical calcu¬ 
lation are listed in Table Q] except otherwise specified. 


A. Shubnikov de Haas oscillation 


where gp r is the Frohlich coupling constant, erfc is the 
complementary error function, and <7 is the effective 
width of the electronic envelope function. 

For MoS 2 on a substrate, the surface optical phonons 
(SOPs) couple to the electrons via an effective electric 
field, which may play important role in transpor t 35 ’ 39 
just like graphene ] 34 ’ 40 ’ 41 The coupling matrix element is 
expressed as^ 


\M TT (q,SO )\ 2 


£ Hr,so f f 
2 £ 0 nq yl + K“ 


1 

1 + 


e ~ 2qd , 


(30) 

with Hr,so the frequency of SOP, and n™ (k°) denoting 
the high (low) frequency dielectric constant of substrate. 

Unlike the case of the static impurity scattering, the 
carrier screening for phonon scattering is dynamic, i.e., it 
is the screening function e(q, n g „) at phonon frequency 


In this subsection we consider the magnetotransport 
of suspended M 0 S 2 at low temperatures. First, the im¬ 
purity scattering is assumed to be only the intravalley 
Coulombic scattering (d = 0). In Fig.[TJ the longitu¬ 
dinal magnetoresistance R xx is calculated versus aver¬ 
age filling factor v 0 = w“ 2 [e| — A 2 ] with e F denoting 
the Fermi energy. The electron density is set to be 
N = 7 x 10 12 cm' 2 and the zero-field linear mobility 
/io = 4000cm 2 /Vs at zero temperature. This value of 
mobility, though one order larger than those currently 
obtained experimentally ) 11 1 13 is consistent with the the¬ 
oretical work— Higher linear mobility at low temper¬ 
ature can be achieved via the gate dielectric engineer¬ 
ing to effectively screen charge impurities,— and doping 
and strain modulations already realized a mobility higher 
than 1000 cm 2 /Vs at room temperature^ we thus expect 
this zero-temperature mobility will be reached in the near 
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TABLE I: Material parameters for monolayer M 0 S 2 used for 
calculation. The T /K subscripts represent intra/intervalley 
phonons. 


Parameter 

Symbol Value 

Lattice constant 

a 

3.193 A 

Hopping integral 

t 

1.1 eV 

Energy gapi 

A 

1.9 eV 

Spin splitting energyi£ 

A 

75 meV 

Mass density^ 

P 

3.1 x 10~ 7 g/cm 2 

Effective Layer thickness* 

a 

4.41 A 

dielectric constant of MoSz^ 
dielectric constant of ZrO^i 

k ; 

7.6 

low frequency 

K° 

24 

high frequency 

K? 

4 

Transverse sound velocity^ 

VTA 

4200 m/s 

Longitudinal sound velocity^ 
Acoustic deformation potentials^ 

ULA 

6700 m/s 

TA 

StA 

1.6 eV 

LA 

Sla 

2.8 eV 

TA 

Dk.ta 

5.9 eV 

LA 

Dk,la 

3.9 eV 

Optical deformation potentials^ 
TO 

D 1 

^r,TO 

4.0 eV 

TO 

Dk,to 

1.9 eV 

LO 

^K,LO 

2.6 x 10 s eV/cm 

Homopolar 

Frohlich coupling^ 

n° 

^r,HP 

4.1 x 10 s eV/cm 

LO 

flFr 

286 meVA 

Phonon energies^ 



TA 

Dk,ta 

23.1 meV 

LA 

Dk,la 

29.1 meV 

TO 

Dr, to 

48.6 meV 

TO 

Dk,to 

46.4 meV 

LO 

Dr,LO 

48.0 meV 

LO 

Dk,lo 

42.2 meV 

Homopolar 

SOP energies of ZrOz^i 

Dr, hp 

D (1 l 

“r.so 

d (2) 

1 J r,so 

50.9 meV 

1st mode 

25.02 meV 

2nd mode 

70.8 meV 


future. 

As can be seen from Fig. [l] the magnetoresistivity ver¬ 
sus filling factor uq or magnetic field B, exhibits marked 
SdHO with a beating pattern, having approximate period 
Av 0 ~ 1 at large fillings or low magnetic fields. The re¬ 
sistivity peaks or valleys locate at integer fillings. There 
is a phase inversion, i.e., a change from the integer fillings 
for peaks to the ones for valleys. These features are in 
vivid contrast to graphene ; 46 i 47 where the SdHO valleys 
locate in the vicinity of half-integer filling factors with¬ 
out beating patterns, but analogous to the behavior of 



FIG. 2: (Color online) Curves of analytical expressions for 
magnetoresistance versus average filling factor vq. The blue 
solid line is obtained from Eq. (El, while the red solid line is 
directly calculated from Eq. (1321) for the maximum fe max = 2. 
The other parameters are the same as Fig.[T| 


conventional 2D electron gas with spin-orbit coupling^ 
Further, at large magnetic fields or small fillings, the pe¬ 
riod of oscillation in monolayer M 0 S 2 halves, which can 
be seen clearly in the inset of Fig.[T| With an increase of 
temperature, the amplitude of SdHO decreases rapidly. 

At these low temperatures shown, contribution to the 
frictional force mainly originates from electron-impurity 
scattering, and thus resistivity R xx ~ —f e i/(N 2 e 2 v). 
For (5-form intravalley short-range scattering potential 
U TT {q) = uo, the zero-temperature magnetoresistivity 
R xx can be expressed as 


_ mul ^ 2 dW 0 Z(q,u) 

Rxx — a rt -) / , 9: 


N 2 e 2 ** Q U 
qr 


uj—0 


N 2 e 2 ll ^ 


gr S {£F)g rs' ( £ f) 


r*+oo 


dzzC ; 


T^t 


(*) 


IDo 


in which the density of states of electrons in the rth 
valley with spin s at Fermi energy ep, 5ts( £ f) = 

_ E n ^ m ^ ns ( £ F)/(27r 2 Zg), can be rewritten, by means 
of Poisson summation formula, as 


5t S ( £ f) = " 2 , 2 { 1 + 2 X! [cos(27rfcz/ TS ) 

tvI-dLu^ 1 


l B w c 


fc=1 


— kp sm(2nki' TS )] exp I —2k 


r 2 


>, (31) 


where v TS = uj~ 2 (e p — A)(ep + A — rsA) is the filling 
factor of electrons in rth valley with spin s, and /3 = 
7 rT 2 /w 2 . For monolayer M 0 S 2 even with low mobility 
/i ~ 10 cm 2 /Vs, the coefficient (3 <C 1, therefore, the 
term with sine function could be omitted safely. In the 
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case of high filling factor v TS , the integral 

r+oo 

/ dzzC. (z) = 2 v Ta S aa >, 

Jo 

and the linear magnetoresistance can be written as 


R 


XX 


riiUQ £p 

2TrN 2 e 2 i£a 4 t 4 ^ TS 

& rs 


^max — OO 

1 + 2 cos(27 rkv TS ) 

fc=l 


x exp 



2 


(32) 


Usually, on the account of the rapid decay of the expo¬ 
nential function, one only needs to keep terms with k = 1 
in the summation, leading to 



_ rijul e|(e| - A 2 ) 
N 2 e 2 na 6 t 6 


1 + 4 cos (2-kvq) 


x cos 



£f +A ) 



(33) 


This represents that the amplitude of oscillation is mod¬ 
ulated by the second cosine function due to the spin¬ 
splitting and there are nodes at A^ 0 /(£f + A) = l ± 
1/4 with l being an integer. Note that three smallest 
nodes in positive regime corresponds to \vq / (e+ + A) = 
0.25,0.75,1.25 or vq = 6.4,19.1,31.9, in agreement with 
the numerical calculation (see Fig-IH). However, the os¬ 
cillating peaks at large magnetic field or small filling fac¬ 
tor obey A(v 0 ) — 0.5 in the figure, which cannot be ex¬ 
plained by the above equation and is due to terms of 
higher frequency. Because of the small value of /3, the 
product (/3 £f/U) 2 may not be considerably larger than 
one and the oscillating terms with k > 1 also may some¬ 
what contribute to the total resistivity. Fig. [2] demon¬ 
strates the results from the approximate expression (l33l) 
and from (1521) with k summing up to 2. It is clear that the 
oscillation part of high frequency comes from the terms 
with k = 2. It is noteworthy that this feature is irre¬ 
spective of the half-integer filling in graphene due to the 
electron-hole symmetry of zero Landau level for massless 
electrons. In the absence of magnetic field, the resistivity 
R xx reduces to 


_ mul £|(£r|-A 2 ) _ mul na 2 t 2 N + A 2 
0 ~ N 2 e 2 it a 6 t e ~ Ne 2 aH 4 ' 1 j 

Despite the linear dispersion on momentum, this resistiv¬ 
ity depends on the electron density owing to its massive 
property, in contrast to the result of graphene.— The 
corresponding density N is 58 x 10 12 cm' 2 when 7 ra 2 t 2 N 
equals A 2 for the present parameters. Hence, for small 
density resistivity Rq is inversely proportional to density 
similar to the case of conventional 2D electron gas, while 
for very large density Rq becomes independent of electron 
density. 


FIG. 3: (Color online) Normalized magnetoresistance versus 
filling factor uo for the case A = 54meV (black solid line) 
and A = 0 (red dash line) at temperature T = 0.3 K. Here 
electron density N = 9.69 x 10 12 cm -2 , linear mobility /ro = 
600 cm 2 /Vs, and or = 8. The inset shows the corresponding 
experimental results, which are replotted as function of vq. 


To compare our theoretical result with recent experi¬ 
mental observation, 2 - Fig. [3] presents the normalized re¬ 
sistivity versus filling factor for another monolayer M 0 S 2 
with electron density N = 9.69 x 10 12 cm -2 and linear 
mobility fi 0 = 600cm 2 /Vs at T = 0.3 K. The curves for 
the cases with spin-orbit coupling A = 54 meV and with¬ 
out spin-orbit coupling are plotted, respectively, as solid 
and dash lines. Here the factor ar = 8 . In the inset, 
we replot the experimental result taken from Fig. 4(b) in 
Ref. ( 23 , as a function of average filling factor. For the ex¬ 
perimental sample having low mobility, the Landau level 
broadening is so large that the beating pattern can not 
be observed. Nevertheless, the phase inversion of SdHO 
peaks still shows clearly. As can be seen, vq = 5 cor¬ 
responds to a position of SdHO peak, while uq = 12 is 
for valley. The numerical calculation agrees with the ex¬ 
perimental observation well. The red dash line for the 
case of A = 0, where the peaks always locate at integer 
filling factors, is also plotted for comparison. We can see 
that the spin-orbit splitting is very important for mag- 
netotransport in monolayer M 0 S 2 , even for low-mobility 
sample in which the full beating pattern of SdHO is not 
easy to observe. 

To investigate the intervalley scattering effect on the 
SdHO, in Fig.[4]we plot the oscillating magnetoresistance 
induced solely by the short-range intervalley disorder at 
various lattice temperatures T = 0.5,1.0,1.5,2.0, 8.0 K. 
Here the relaxation time r s = l/(m*riiufy is set to be 
1 ps. For the purpose of comparison, the SdHO induced 
by intravalley short-range electron-impurity scattering is 
also plotted in the inset of Fig. 01(a) for the same value of 
relaxation time at T = 0.5 K. It is seen that the magne- 




























































FIG. 4: (Color online) (a) Intervalley electron-impurity scat¬ 
tering induced magnetoresistance vs the average filling factor 
vo at different lattice temperatures T = 0.5,1.0,1.5, 2.0, 8.0 K 
when the relaxation time t s = 1 ps. The inset shows the 
intravalley short-range electron-impurity scattering induced 
magnetoresistance at T = 0.5 K for the same relaxation time, 
(b) Magnetoresistance contributions , R^ x and R^J 

versus the average filling factor at T = 0.5 K. 


toresistance induced by the intervalley collision, though 
almost two order smaller than intravalley one, also ex¬ 
hibits SdHO versus the average filling factor and the ex¬ 
trema show up at integer fillings and the oscillation is 
also modulated by the spin-orbit interaction with nodes 
locating at the same positions as in the intrasuband case. 
But the modulation appears much stronger than the in¬ 
travalley one: with increasing temperature the amplitude 
of SdHO decreases, while the envelope of oscillation still 
exists even at T = 8.0 K when the intravalley one dis¬ 
appears. Especially, in contrast to the intravalley case, 
there exists a regime AB (3 < v 0 < 8 or 9 T < B < 24 T), 
in which the magnetoresistance almost vanishes. 

All these can be referred to the fact that, in contrast 
to intravalley case, the intervalley scattering hardly takes 
place between two states having the same Landau in¬ 
dex n > 0. The Landau levels e Tns expressed in © for 
n = 1,2,3,..., can be written as with i = ts. As 
indicated in Eq. (1211) the resistivity is proportional to the 
product of DOSs of two close (contributory) Landau lev¬ 
els around the Fermi energy. In the vicinity of Fermi 


energy, for a fixed Landau index n the level separation 
of different l is almost independent of the magnetic field, 
while the distance between Landau levels having same 
l but different Landau indexes n and n' is proportional 
to the magnetic field. Hence, at large magnetic fields 
two contributory Landau levels of different l must have 
the same Landau index n. At low magnetic fields, the 
Landau indexes of two contributory Landau levels may 
not be equal to each other and their difference increases 
with decreasing magnetic field. In Fig.QJb) the magne¬ 
toresistance R^\ contributed from electron transitions 
between two Landau levels with Landau-index difference 
of to near Fermi energy, are plotted as functions of av¬ 
erage filling factor vq at 0.5 K. At low filling factors or 
large magnetic fields, the energy distance between levels 
with same Landau and spin indexes but different valley 
indexes is smaller compared with that between Landau 
levels with different Landau indexes, and we only need to 
consider the transition between levels of different valleys 
but having same Landau index vq, leading to R X J ■ For 
low magnetic fields, contributions of electron transitions 
between levels having different Landau indexes dominate. 
Here, R,i] r - stands for contribution from the transitions 
between vq and vq — 1 levels and those between vq and 
vq + 1 levels. Rxx stands for contribution from electron 
transitions between vq + 1 and vq — 1 levels, and Rxx for 
contribution from transitions between vq + 2 and vq — 1 
levels and those between vq — 2 and vq + 1 levels. It is 
found that R^ x + R^ x + R'xx + R*xx almost equals the 
total magnetoresistance R xx shown in Fig. [2(a). 

Rxx becomes quite small when e„ 0j+ — ep > F and/or 
£f — e "o,- > r, i-e., it is almost zero for magnetic fields 
lower than a certain value. On the other hand, with the 
increase of the magnetic field, the level distance of dif¬ 
ferent Landau indexes enlarges, leading to R^ x almost 
vanishing for magnetic fields larger than a certain value, 
which is determined by — £ uo -i,+ > L (so do for 
RxJ and r£ 3 J). In the range between these two magnetic 
fields, the total magnetoresistance appears very small. 
For the present parameters (set e^ 0l + — £f = 1.3r), this 
range is 8.8 T < B < 23.3 T or 3.1 < vq < 8.2, as indi¬ 
cated AB in the Fig.[I] 


B. Magnetophonon resonance 

Now we concentrate on the case of higher tempera¬ 
ture up to room temperature. First, we consider the 
suspended M 0 S 2 . The total magnetoresistances R xx in¬ 
duced by the intravalley screened Coulombic electron- 
impurity scattering (d = 0) and all above-mentioned in¬ 
travalley and intervalley electron-phonon couplings ex¬ 
cept SOP mode, are plotted as functions of magnetic 
field at various high temperatures in Fig.[5](a). The R xx 
increases with the increment of magnetic field, accom¬ 
panying an oscillation at large fields. The behavior of 
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FIG. 5: (Color online) (a) Total longitudinal magnetoresis¬ 
tance R xx in suspended M 0 S 2 versus magnetic field B at high 
temperatures T = 150, 200, 250, 300 K. The inset shows the 
main contributions from electron-phonon scattering, where 
R oa = R&c + Rop■ (b) The impurity-induced magnetoresis¬ 
tance Rim is plotted as a function of magnetic field B. (c) The 
magnetoresistance induced by intravalley transverse (longitu¬ 
dinal) acoustic phonons Rta(Rla) versus the ratio u 0 /ub. 
(d) The magnetoresistance induced by intravalley homopolar 
(Frohlich coupling) optical phonons Rhp (Rfr) versus the ra¬ 
tio w 0 /cjb at temperatures T = 100, 200, 300 K. The other 
parameters are the same as Fig.|T] 


resistivity increase with increasing magnetic field is due 
to impurity-induced resistivity Ri m as shown in Fig.(5](b). 
Since SdHO almost disappears at this temperature, the 
small oscillation in R xx originates from phonon scatter¬ 
ings. With ascending temperature, R lm descends, while 
the total R xx increases because of the increasing con¬ 
tributions from electron-phonon scatterings. It is found 
that, in addition to the electron-impurity scattering, the 
contributions of the intravalley transverse and longitu¬ 
dinal acoustic phonons, Rta and Rla, and those of ho¬ 
mopolar and Frohlich coupling optical phonons, Rhp and 
Rfr, play a dominant role in the total resistivity. The 
inset of Fig. [5j( a) shows that the oscillation arises mainly 
from the optical contribution R op = Rhp + Rfr- The 
acoustic one R ac = Rta + Rla gives almost a constant 
value at room temperature. 

Actually, the acoustic contributions Rta, Rla also os¬ 
cillate with magnetic field, especially at relatively low 
temperature, exhibiting the so called MPR induced by 
acoustic phonons , 34 ' 51- — The magnetoresistance peaks 


occurs when the energy of the optimum phonons w a = 
2 &FU ac equals an integral multiple of the inter-Landau- 
level distance ub near Fermi surface. Here v ac = vta 
or ula is the sound velocity for the transverse or lon¬ 
gitudinal mode. The energy distance ub between two 
intravalley Landau levels with same spin around Fermi 
energy £p is given by 

2 2 

wb ~ — , ~ — • (35) 

2 \J (A - TS \) 2 + l^rs^c ^ 

Fig-EKc) indeed shows the oscillation of magnetoresis¬ 
tance for both Rta and -Rla with inverse magnetic field 
having period A(u> a /u>B) — 1. With increasing temper¬ 
ature, the peaks at high ratio w 0 /ub tend to disappear 
gradually. Further, the peak slightly shifts to smaller 
uia/ijJ b position and the magnetoresistance due to longi¬ 
tudinal mode Rla becomes larger than the transverse one 
Rta in view of enlarged phonon energy with the rise of 
temperature. In the present case with relatively low mo¬ 
bility, the MPR induced by acoustic phonons has little 
influence on the total magnetoresistance. However, for 
monolayer M 0 S 2 having ultrahigh mobility, the acoustic 
electron-phonon coupling contributes dominantly at low 
temperature, hence this MPR could be observable. 

With further increase in lattice temperature, the 
electron-optic phonon coupling becomes more and more 
important in comparison with other scattering mecha¬ 
nisms. Due to the large coupling coefficients, the re¬ 
sistivities induced by the homopolar and Frohlich in¬ 
teractions have largest values. It is well known that 
the resistivity exhibits MPR when the energy of opti¬ 
cal phonons uj 0 equals the distance of Landau levels. In 
monolayer M 0 S 2 for £p — A <C A, the intravalley Landau 
levels are almost evenly spaced and the level distance 
approximately equals ub- In Fig.[5](d) the magnetore¬ 
sistance Rrp or Rfr is plotted versus the ratio w 0 /ub, 
where lo 0 = Df,hp or Df,lo is respectively the frequency 
for homopolar or Frohlich coupling. It is true that the 
magnetoresistances show peaks or valleys at w d /ub = l, 
i.e. magnetoresistance oscillates with inverse magnetic 
field having period A(uj 0 /lob ) — 1. However, in con¬ 
trast to the usual MPR induced by optical phonons in 
two-dimensional electron gas, the oscillating resistivity 
in M 0 S 2 is modulated due to the spin splitting by an 
approximate factor cos ^27 t^- analogous to the 

SdHO. Hence, there are nodes at —— 7 -r = l ± 7 . This 
leads to the nodes appearing at w„/ub = 19.1,31.9, in 
accordance with Fig.[5](d). 

Now we study the MPR for M 0 S 2 on a Zr 02 sub¬ 
strate. It is found that the frequencies of SOPs for 
Zr 02 are so small that they play an important role in 
electron transport . 3 ' 5 Hence, in Fig. [G] magnetoresistances 
for M 0 S 2 on Zr 02 are plotted versus magnetic held. In 
the calculation, the elastic scattering is assumed to be 
the intravalley remote-impurity scattering distributing at 
d = 1 nm from the single layer, the inelastic scatterings 
are due to all the intrinsic modes mentioned above and 
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FIG. 6: (Color online) Magnetoresistance in ZrC> 2 /MoS 2 /Air 
structure at high temperature T = 200, 300 K. Here the elec¬ 
tron density N = 7 x 10 12 cm -2 . In (a) and (b), the zero- 
field mobility at zero temperature fj ,o = 4000 cm 2 /Vs and 
or = 3, while in (c) and (d) fj,o = 600 cm 2 /Vs and ar = 8 
for another sample, (a) Total magnetoresistance R xx and 
impurity-induced one versus magnetic field B. (b) The 

SOP-induced magnetoresistance Rqq is plotted as a function 
of fipjo/wB- (c) The total magnetoresistance for the sample 
with low mobility versus the magnetic field, (d) The magne¬ 
toresistance induced by intravalley homopolar, Frohlich cou¬ 
pling optical phonons, and SOP Rhp,Rfk, and R^q versus 
magnetic field. 


the intravalley SOPs. (/j,® used in the screening of elas¬ 
tic scattering is estimated to be 0.3<7tf from Fig. 2 in 
Ref.[H. Two samples with different zero-field mobilities 
are considered for comparison. For the clean system, the 
remote impurity scattering weakens the quick increase of 
impurity-induced resistivity with magnetic field in con¬ 
trast to suspended case, leading to more evident MPR in 
the total magnetoresistance. On the other hand, in com¬ 
parison to suspended M 0 S 2 , the MPR behavior becomes 
more complex because of the crucial influence of SOPs, 
especially the mode with low frequency f2p^ 0 . The resis¬ 
tivity Rqq induced by the first SOP mode is plotted in 

Fig-Etb) versus H^ so /wb. The resonant feature of i?gg 
is similar to other intrinsic modes. However, for another 
sample with low mobility in heavily overlapping-Landau- 
level regime, the MPR almost disappear. The magne- 
toresistance increases monotonously with magnetic field, 


and only a small oscillation occurs at very large field. In 
Fig. [10(d), contributions of three important optical modes 
are plotted. 

The effect of a SiC >2 substrate on the MPR of M 0 S 2 
is also tested and it is found that this dielectric plays 
negligible role due to its large frequencies of SOPs. 


IV. SUMMARY 

In summary, we have studied the linear magnetotrans¬ 
port in single layer M 0 S 2 employing a balance equa¬ 
tion analysis by including spin-orbit coupling and all 
kinds of intravalley and intervalley electron-impurity and 
electron-phonon scatterings. 

The existence of an energy gap between the conduction 
and valence bands, or lack of electron-hole symmetry of 
the zero Landau level in M 0 S 2 , makes its magnetotrans¬ 
port behavior more like a conventional 2D electron gas 
than graphene: the resistivity peaks or valleys of its low- 
temperature SdHO, resulting either from intravalley or 
from intervally clastic scatterings, locate at integers of 
filling factor vq. 

The large spin-orbit coupling in the system, however, 
gives rise to a significant modulation or beating of the 
magnetoresistance oscillation, or a phase inversion of the 
oscillation peaks. The agreement between theoretical 
prediction and recent experiment on the phase inversion 
of SdHO peaks demonstrates the importance of the spin- 
orbit splitting in magnetotransport even for systems of 
low-mobility. The clear beating pattern of the oscillating 
magnetoresistance should appear in the well-separated 
Landau-level regime in high-mobility systems. 

On the other hand, the behavior of magnetoresistance 
oscillation at large magnetic fields or small filling factors 
appears different for intravalley and intervalley scatter¬ 
ings: the period of oscillation associated with intraval¬ 
ley scattering may halve due to the weak decay of the 
second-order oscillating term, while in the case of inter- 
valley disorder much stronger spin-orbit induced SdHO 
modulation shows up that there exists a magnetic-field 
range in which the magnetoresistivity almost vanishes. 
Of course, intervalley elastic scattering contributes only 
a much smaller part to the total magnetoresistance than 
that from intravalley ones. 

At high temperatures, the magnetoresistance oscilla¬ 
tion arising from MPR may show up in the smooth 
impurity-induced resistivity background both for sus¬ 
pended and nonsuspended samples with high mobility. 
Both acoustic phonons (mainly intravalley transverse and 
longitudinal acoustic modes) and optic phonons (mainly 
homopolar and Frohlich modes) can induce MPR. A 
beating pattern with the same frequency as in the SdHO 
also appears in the optical-phonon-induced MPR due to 
spin-orbit coupling. For the single layer on a substrate, 
another resonance due to SOPs may occur, resulting in a 
complex behavior of the total magnetoresistance. How¬ 
ever, for nonsuspended layer with low mobility, the mag- 
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netoresistance oscillation almost disappears and the re¬ 
sistivity increases with field monotonously. 
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